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The small-angle scattering curves of deterministic mass fractals are studied and analyzed in the momentum 
space. In the fractal region, the curve I(q)q D is found to be log-periodic with a good accuracy, and the period 
is equal to the scaling factor of the fractal. Here D and I(q) are the fractal dimension and the scattering 
intensity, respectively. The number of periods of this curve coincides with the number of fractal iterations. We 
show that the log-periodicity of I(q)q D in the momentum space is related to the log-periodicity of the quantity 
g(r)r 3 ~ D in the real space, where g(r) is the pair distribution function. The minima and maxima positions 
of the scattering intensity are estimated explicitly by relating them to the pair distance distribution in the real 
space. It is shown that the minima and maxima are damped with increasing polydispersity of the fractal sets; 
however, they remain quite pronounced even at sufficiently large values of polydispersity. A generalized self- 
similar Vicsek fractal with controllable fractal dimension is introduced, and its scattering properties are studied 
to illustrate the above findings. In contrast with the usual methods, the present analysis allows us to obtain not 
only the fractal dimension and the edges of the fractal region, but also the fractal iteration number, the scaling 
factor, and the number of structural units from which the fractal is composed. 
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I. INTRODUCTION 

Modern experimental techniques in materials science for 
preparing hierarchically organized systems at nano and micro 
scales open up new possibilities in controlling their functions 
and properties. Considerable scientific and technological ef- 
forts have been directed to the development of such systems 
ll Ullll . The hierarchical structures, like fractals (12=] or bio- 
logical objects HH], manifest themselves in the electromag- 
netic 0], statistical 0], dynamical d] or optical fnlfiai 
properties, and therefore one of the basic challenges is to un- 
derstand correlations between these properties and material's 
microstructure 119ll20il . 

Experimentally, the material's microstructure can be deter- 
mined by means of scattering techniques, which yield elastic 
cross section per unit solid angle as a function of momentum 
transfer. The cross section, usually normalized per unit vol- 
ume of a sample, is called scattering intensity: I(q) = . 

V ail 

The scattering angle 2-d is related to the momentum transfer 
Hq by equation q = AttX^ 1 sini?, where A is the radiation 
wavelength. The typical range in which the materials exhibit 
a hierarchical structure is within d = 1 -j- 1000 nm (see, 
e.g., Ref. Il2ll0 . The structural characteristics in this range 
are probed well when the scattering wave vector lies within 
q ~ 2n/d= 10~ 3 -^1 A , which corresponds to small-angle 
X-ray or neutron scattering (SAXS/SANS) lHHH. This im- 
portant investigation technique has become a powerful tool for 
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studying the fractal microstructure (see reviews 1E5U28I1 and 
references therein). 

Presently, various deterministic fractals can be artificially 
created Jj-Ja, [H [l(J due to a rapid progress in nanotech- 
nologies. This makes it relevant and promising to study the 
connections between their scattering properties and their mi- 
crostructure. Deterministic fractals frequently allow an an- 
alytical description of small-angle scattering (SAS) and thus 
give us "exactly solvable models", very useful for understand- 
ing the scattering properties of fractals. By choosing parame- 
ters of deterministic fractals at random (say, by introducing 
polydispersity), one can understand the basic properties of 
random fractals as well [29, 30]. 

An essential feature of fractal is Hausdorff (fractal) dimen- 
sion. While it can be rigorously defined (see AppendixlAl, it 
is convenient in practice to ado pt a simple descriptive defini- 
tion of the dimension D HHHHH: N oc {l/a) D for a ->• 0, 
where N is the minimum number of open sets of diameter a 
needed to cover a fractal. Non-random (deterministic) frac- 
tals are generated by deterministic processes, in particular, by 
iterative rules. This implies the presence of an initial set (ini- 
tiator) and a generator (iterative operation). The number of 
iterative operations is called fractal iteration. Let us consider 
a finite iteration of mass fractal, consisting of simple units, 
say, balls of radius a. Let the fractal size be I. If N is the total 
number of balls then the fractal dimension D is defined by the 
asymptotics 

Ncx(l/a) D (1) 

for a large number of iterations, which assumes a — > (from 
the physical point of view, a <C 0- One can use this defi- 
nition for a given finite iteration (a = const) in order to es- 
timate the number of balls enclosed by an imaginary sphere 
of radius r with a ball in the center. Equation ([TJi yields 
N(r) oc (r/a) D oc r D . This estimation is valid within the 
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fractal region l mm < r < I, in which fractal properties can be 
observed experimentally. Here Z m ; u is the minimal distance 
between the ball centers. 

A main indicator of the fractal structure is the power-law 
dependence between the small-angle scatte ring intensity and 
the absolute value of scattering vector | t 34 - 13711 . For a mass- 
fractal, it takes the form 

I(q) a q-°. (2) 

The fractal region in the real space implies the fractal region 
in the reciprocal space 

1A < 9 < 1/imin, 0) 

for which Eq. © is applicable. Equation (O is quite general, 
because the fractal dimension essentially regulates the spa- 
tial correlations between fractal units. Indeed, the scattering 
intensity is proportional to the Fourier transform of the pair 
distribution function g(r), describing the spatial correlations 
between particles inside a fractal (see Sec. I VI Al below). In 
more detail, once a particle lies in the coordinate origin then 
the number of other particles in the volume d 3 r near the radius 
vector r is equal to dN = ng(r)d 3 r, where n is the average 
particle density. If the pair distribution function is radially 
symmetric, we obtain N(r) = J Q ng(r')4:Trr' 2 dr' for the to- 
tal number of particles in the sphere of radius r. On the other 
hand, the fractal dimension implies N(r) oc r D , as discussed 
above. It follows that 

g(r) cc r°- 3 (4) 

for r < I lf38i[39tl . which leads to Eq. (0 in accordance with 
Erdelyi's theorem for asymptotic expansion of Fourier inte- 
grals Ml . 

Thus, three parameters can be extracted from experimen- 
tal fractal scattering intensities: the exponent and the edges of 
the fractal region in the q-space, which appear as "knees" in 
the scattering line on a logarithmic scale. Other parameters 
are not usually obtained from SAS curves. For random frac- 
tals, it is rather difficult to extract more information from the 
scattering, because a fine structure of particle correlations is 
usually smeared due to the randomness. One can expect that 
deterministic fractals, being more ordered, allow us to obtain 
additional information from the scattering data. In this paper 
we suggest a scheme for estimating the iteration number, the 
scaling factor and the total number of structural units in deter- 
ministic mass fractals from the scattering curves. 

The present analysis of deterministic fractals is based on 
the method of calculating the scattering amplitude, devel- 
oped in the papers |Jj| [30, SH l42ll (for the generalization, 
see Ref. B43IP . Recently, analytical calculations of SAS in- 
tensity from the 3D triadic Cantor and Vics ek II44I1 sets were 
reported in Ref. j45tl . In the previous paper ll46ll . the general- 
ized three-dimensional (3D) Cantor set was suggested, whose 
dimension can vary from to 3 by means of the scaling fac- 
tor. If the SAS is considered from monodispersive sets, which 
are randomly oriented and placed, then the scattering intensi- 
ties represent minima and maxima superimposed on a power- 
law decay, with the exponent equal to the fractal dimension 



of mass fractal. This dependence of the intensity on momen- 
tum is called the generalized power-law decay. As was shown 
l42l l45l l46ll . the minima and maxima are damped with in- 
creasing polydispersity of the fractal sets. The physical rea- 
sons for such behavior are quite clear: the fractal dimension 
dictates the power law (|4]i for the particle correlations only on 
the average. In this paper, we show that in a deterministic 
fractal, the pair distribution function g(r) also obeys the gen- 
eralized power law [see Eq. d55l ) below], that is, it exhibits 
minima and maxima on the power-law decay ©. This struc- 
ture appears due to clusterization and intimately relates to the 
fractal scaling factor. Polydispersity smears the spatial dis- 
tribution between fractal units. Strongly developed polydis- 
persity thus leads to the simple power-law behavior and 
©. In this paper, we construct and consider the generalized 
self-similar Vicsek fractals (GSSVF) as an example in order 
to elucidate the above scattering properties. 

The paper is organized as follows: in Sec. [II] we emphasize 
some important issues concerning SAS. Section Hill describes 
the construction of GSSVF with controllable dimension, gov- 
erned by the scaling factor. In the subsequent section, we 
derive analytically the scattering amplitude for GSSVF and 
calculate its scattering properties: the intensity, structure fac- 
tor, and radius of gyration. The influence of polydispersity on 
the fractal scattering properties is considered in Sec. [V] The 
spatial correlation functions of the fractal are obtained and in- 
terpreted in Sec. [VI] In Conclusion we discuss the obtained 
results and promising prospects of the developed analysis. 



II. GENERAL REMARKS ON SMALL-ANGLE 
SCATTERING 

Let us consider the SAS scattering (neutron, X-ray, light, 
or electron diffraction) on a sample consisting of microscopic 
objects with the scattering length bj. If the multiscattering 
processes are neglected (which is a very good approxima- 
tion usually), then the differential cross section of the sam- 
ple is given by JH] dcr/dfi = |^4(<z)| 2 , where A(q) = 
J V t p s (r)e tqr d 3 r is the total scattering amplitude, and V' is 
the total volume irradiated by the incident beam. The scat- 
tering length density can be defined with the help of Dirac's 
(^-function as p s (r) — J2j bjS(r — rj), where Tj are the mi- 
croscopic object positions. 

In this paper, we restrict ourselves to two-phase systems, 
which are composed of homogeneous units of "mass" density 
p m . The units are immersed into a solid matrix of "pore" den- 
sity pp. A constant shift of the scattering length density in the 
overall sample is important only for small values of wave vec- 
tor q < 2-7T / (V 1 ) 1 / 3 , which are usually beyond the resolution 
of scattering device. Therefore, by subtracting the "pore" den- 
sity, we can consider the system as if the units were "frozen" 
in a vacuum and had the density Ap = p m — p p . The density 
Ap is called scattering contrast. 

In practice, it is convenient to represent the total scattering 
amplitude as a sum of amplitudes of rigid objects. For in- 
stance, considering the scattering from stiff fractals, whose 
spatial positions and orientations are uncorrected, one can 



choose them as the objects. Then the scattering intensity (that 
is, the differential cross section per unit volume of the sample) 
is given by 



I(q) = n\Ap\ 2 V 2 ( \F(q 



(5) 



where n is the fractal concentration, V is the volume of each 
fractal, and F(q) is the normalized formfactor 



(6) 



obeying the condition F(0) = 1. Here, the brackets (• • • ) 
stand for the ensemble averaging over all orientations of the 
fractals. If the probability of any orientation is the same, then 
it can be calculated by averaging over all directions n of the 
momentum transfer q = qn, that is, by integrating over the 
solid angle in the spherical coordinates q x = q cos (p sin -d, 
q y = q sin tp sin it) and q z = q cos d 
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(f(qx,q v ,qz)) = J dtisiwd J dp> f(q,$,<p). (7) 
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The intensity in zero angle results directly from Eq. (0 

1(0) = n\Ap\ 2 V 2 . 



(8) 



Given the contrast, concentration of fractals, and the absolute 
value of intensity, the fractal volume can be determined from 
Eq. ®. 

Once a deterministic fractal is composed of the same ob- 
jects, say, N balls of the same radius R, then the formfactor 
can be written as 

Hi) = P q F (qR)/N. (9) 

Here the ball formfactor of unit radius is given by ll23ll 

F (z) = 3(sin z — z cos z) / z 3 , (10) 

and p q = J^j e~ lq rj is the Fourier component of the density 
of ball centers, where rj are the center-of-mass positions of 
balls. By substituting Eq. (0 into Eq. (O, we obtain for the 
scattering intensity 



I(q) = I(Q)S(q)\F (qR)\ 2 /N, 
where we put by definition (see, e.g., l47ll ) 
S(q) EE ( PqP - q )/N. 



(11) 



(12) 



The quantity S(q) is called structure factor, and it is intimately 
connected to the pair distribution function (see the discussion 
in Sec. [VI] below). The structure factor carries information 
about the relative positions of the balls in the fractal and can 
be rewritten as 



1 N 



(13) 






FIG. 1. (Color online) The initiator and first three iterations for the 
GSSVF at the scaling factor /3 B = 1/6. (See Ref. |48] for a video of 
the first iteration of GSSVF.) 



It obeys the relations 5(0) = N and S(q) ~ 1 for q > l// m i„. 
The former relation follows from the definition ( TT2b . and the 
latter is fulfilled, because the contribution of non-diagonal 



terms in the r.h.s. of Eq. (O Y,jjtk {e~ iq ' {r ' ~ Tk ' ) /N tends to 
zero in the limit of large momentum due to the randomness of 
the phase. This limit is analogous to the limit of geometrical 
optics. 

Let us emphasize that although Eq. (TTTT) is quite general, the 
choice of structure factor S(q) and the formactor Fq is rather 
arbitrary. It depends on the choice of the scattering units, to 
which the formfactor Fq is related. This choice means just 
a regrouping of terms in the total scattering amplitude A(q), 
introduced in the beginning of this section. For example, if 
we associate Fo with the formfactor of entire fractal, then the 
structure factor describes the spatial correlations between dif- 
ferent fractals. In this case, the structure factor equals one for 
non-zero momenta, because the fractal positions are assumed 
to be completely uncorrelated. For this choice, Eq. © is anal- 
ogous to Eq. ( TTTT l with S(q) = 1. 

An important characteristic in SAS is radius of gyration i? g . 
It can be defined from the intensity expansion in the Guinier 
regime q < l/l B23I1 



I(q)=I(Q)(l-q 2 R 2 J3 + •■■)• 



(14) 



III. 



CONSTRUCTION OF GENERALIZED SELF-SIMILAR 
VICSEK FRACTAL 



The construction of the generalized self-similar Vicsek 
fractal (GSSVF), embedded into three-dimensional space, is 



4 



very similar to that of generalized Cantor fractals considered 
inRef. JH. 

We start with a cube with the edge Iq and place in its center 
a ball of radius Iq/2. This is zero fractal iteration, called the 
initiator. Let us choose the Cartesian coordinates: the origin 
lies in the cube center, and the axes are parallel to the cube 
edges. The iteration rule (generator) is to replace the initial 
ball by nine smaller balls of radius (3Jq/2. The position of 
one ball is at the origin, while the centers of the eight other 
balls are shifted from the origin by the vectors 



aj ={±(3 t l ,±p t l ,±p t l } 
with all the combinations of the signs, where 

& = (1 - ft)/2. 



(15) 



(16) 



The dimensionless positive parameter (3 S , called scaling fac- 
tor, obeys the condition j3 s < y3/ (v3 + 2) in order to avoid 
overlapping of the balls. Next iterations are obtained by per- 
forming an analogous operation to to each ball of the first it- 
eration, and so on (see Fig. Q]). Infinite number of iterations 
yields the ideal GSSVF. 

Thus, the total number of balls at the mth iteration equals 



and the corresponding radius is given by 



& m V2. 



(17) 



(18) 



In accordance with the definition ((TJ, the fractal dimension 
can be calculated as (see also AppendixlAt 



D = lim 



loeiV„ 



m^oo log(7 /r„ 



2 log 3 
'logA' 



(19) 



Because of the restrictions imposed on the scaling factor, the 
fractal dimension of GSSVF can vary from to 2.862 . . .. 

The above procedure is an obvious generalization of the 
original Vicsek fractal B44I1 . which is constructed from the 
cubes by means of the scaling factor (3 S = 1/3 and whose 
fractal dimension is equal to 2. 



IV. THE FRACTAL FORM FACTOR AND STRUCTURE 
FACTOR 

A. The analytical expressions 

The standard method for calculating the fractal formfactor 
is the Debye formula j22il23il49tl . including the double inte- 
gration over the scattering length density. Direct application 
of the formula could be extremely difficult, since the number 
of balls in the fractal increases exponentially with the itera- 
tion number. Here we apply an analytical method of calculat- 
ing the fractal formfactor suggested in Refs. 

SUSS! (see 

also Ref. |[50ll . where the Fourier transform of the Cantor mea- 
sure was calculated). This scheme was developed afterwards 
in Refs. H@ |43L l46ll . We give only the results here and 
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FIG. 2. (Color online) (a) Scattering intensity J23t for the first four 
iterations of monodisperse GSSVF. The scattering curve for the mth 
iteration is scaled up for clarity by the factor I0 3 ' m_1 ' . Black at low 
q shows the Guinier regions; light gray (orange) represents fractal 
regions; dark gray (blue) at high q Porod regions, (b) The fractal 
structure factor l |25l > for the first four iterations. Inset: the quan- 
tity S(q){ql Q ) D /N m versus qlo, log-periodic in the fractal region, 
(c) The structure factor l |25l > (full lines) and its approximation J28t 
(dashed lines) for three iterations. The values of structure factor for 
m — 3 and m = 2 are scaled-up by 10 and 100, respectively. 
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refer the reader to our previous paper l46ll for details of the 
derivation. 

Once the number of balls at the mth iteration of GSSVF 
( fTTI i and their radius ( fT8l) are known, the total fractal volume 
is given by 



V m = N m /3f l V , 



(20) 



where Vb = -j-(y) 3 is the volume of initial ball. Then the 
scattering intensity in zero angle I m (0) can be calculated with 
the help of Eq. ©. 

The formfactor of the mth generation is calculated analyti- 
cally by means of the generative function, which is determined 
by the positions of the centers of balls inside the fractal for 
each iteration. For the GSSVF, the generative function reads 

G m (q) = [l + 8cos(u m q x )cos(u m q y )cos{u m q z )]/9, (21) 

with u m = /o/^t/?™ -1 and m = 1, 2, • • • . Here f3 t is defined 
by Eq. ( [TBI . If we put Go(q) = 1, one can express the form 
factor as 



F m {q) = F (r m g) J! Gi(9) 



(22) 



for m = 0, 1, 2, • • • . The form factor of ball and its radius for 
the mth iteration are given by Eqs. ( fTOb and ( fl"8T >, respectively. 
Then SAS intensity is obtained with the help of Eq. (0 by 
averaging over the solid angle © 



I m (q)/I m (0) = (\F m (q)\ 



(23) 

Equation (l23l l gives the normalized scattering intensity as a 
function of qlo and the scaling factor /3 S . 

Comparing Eqs. (O and (l22l . one can conclude that the 
Fourier component of the density of ball centers in the mth 
iteration is given by 

m 

P^=N m llGiiq). (24) 

i=0 

By substituting this expression into Eq. (V2\ , we obtain 

m 

S m (q)/N m = (l[\G i (q)\ 2 y (25) 

It follows from Eqs. ([22}, ([23), and <|25} that 

I m (q)/I m (0) = \F (/3™ql /2)\ 2 S m (q)/N m , (26) 
which is in agreement with the general relation (fTTI i. 

B. Analysis and interpretation of the obtained results 

First, let us discuss the physical meaning of the scattering 
amplitude ( f22b . It can be written in different forms 



= G 1 (q)G 1 (/3 s < 7 )^ m _ 2 (/3 2 < ? ) 



The resulting amplitude can be understood as a sum of am- 
plitudes of different clusters in the fractal, see Fig. [JJ As 
discussed in Sec. HU the choice of clusters is quite arbitrary. 
Thus, the amplitude d2Tb is a sum of N\ = 9 amplitudes 
of (m — l)th fractal iterations, with being the density 
Fourier component of their centers-of-masses. Or it is a 
sum of N2 = 9 2 scattering amplitudes of (m — 2)th itera- 
tions, whose spatial positions are described by p q , and so 
on. Then we can represent the resulting intensity for each 
choice of clusters as a product of associated structure factor 
and formfactor (\F n (q)\ 2 ) ~ £ 1 (g)<|F TO _ 1 (&«i)| a )/JVi ~ 
Si(9)Si(A?)(|f m -i(i3 B 2 ?)P)M 2 ~ ••• ~ S 1 (q)S 1 (f3 s q) 
■ ■ ■ 5i(/3™- 1 q)| J F 1 (/3 s m ^o/2)| 2 /A r m. Therefore, we arrive at 
the approximation for the fractal structure factor 



5 m (g)~5 1 (g)5 1 (/3 s g)---5 1 (/3 s m " 1 '7)- 



(28) 



The numerical results for the scattering intensities and re- 
lated structure factor are shown in Fig. [2] One can see a com- 
plex pattern of maxima and minima superimposed on a power- 
law decay (O, that is, the generalized power-law behavior of 
intensity, discussed in detail below. In the curves, one can 
separate out three different regions: the Guinier, fractal, and 
Porod regions. 



1. Guinier region 

The Guinier region, as discussed in Sec. [Ill is 

q < l/lo- (29) 

In this region, the intensity is well described by the expansion 
(TT4b . The intensity at zero angle is given by Eq. ^ at the 
fractal volume d20i i. and S m (0) — N m . Expanding the form- 
factor (|22"1) in power series in qlo and substituting the result 
into Eq. (l23l l yield the fractal radius of gyration for the mth 
generation 



Ra 



1 ~ 4 2 ' 



-I 2 



(30) 



Here, for a uniform ball of radius Iq/2, the radius of gyration 
is given by R g0 = (3/5) 1 / 2 / /2. In the limit of high number 
of iterations (the ideal fractal), when f3 2,n <C 1, we obtain 



Re 



(8/3) 1 / 2 /3 t / 

(1-#)V2- 



(31) 



2. Fractal region 

The fractal region in the momentum space is determined 
by the maximal and minimal distances between the ball cen- 
ters, see Eq. (0. As can be seen from the fractal construction 
(Sec. HIIli. these distances are of order Iq and Zo/3t/5" l_1 , re- 
spectively. This yields for the fractal region 



(27) 



1/J0 < q < v(w; 4 )' 



(32) 
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The explicit analytical expressions (f2TT>-(f23T> allows us to 
check up these general estimations. For a given momentum, 
if the cosine argument in G m +i is much smaller than 1, then 
G TO +i ~ 1, and further increasing of m does not lead to an 
essential correction. Hence, the mth iteration reproduces the 
intensity that the ideal fractal would give at this momentum. 
As a result, the scattering amplitude of the mth iteration coin- 
cides with the amplitude of the ideal fractal within the region 
i$ 1> which is consistent with Eq. d32b . Such a 
behaviour can be seen in Fig. [2] as the coincidence of various 
iterations in the fractal regions. In the fractal region, the value 
of normalized intensity is very close to that of structure factor 

I m (q)/Im(P) - S m {q)/N m , (33) 

since Fq ~ 1 in Eq. ( f26b . 

One can interpret the structure factor S\ (q) by analogy with 
optics. The quantity NiGi(q) is nothing else but the ampli- 
tude, produced by nine interfering point-like particles of unit 
amplitude. The resulting scattering pattern iVf (|Gi(q)| 2 ) = 
Si(q)Ni is the intensity, averaged over all directions of vec- 
tor q, see the representation ( fT3l . In optics, this average cor- 
responds to diffraction with an entirely uncollimated beam, 
which leads to the strong spatial incoherence. As a result, 
only the first minimum and maximum are quite distinguish- 
able. They are associated with out-of-phase (tp = tt) and 
in-phase (tp = 2ir) interferences, respectively. The next min- 
ima (tp = 3-7T, 57T, • • • ) and maxima (tp = 4ir, 6tt, ■ ■ ■ ) are not 
practically seen; however, they are responsible for the subse- 
quent series of small oscillations. Besides, at q — we have 
a completely coherent diffraction (the intensity is equal to the 
squared number of the points), and at large q, completely in- 
coherent regime of geometrical optics (the intensity is equal to 
the number of the points). By Eq. (13[ , the phase is governed 
by distances between the points, which are of order s/Sptlo, 
and tp ~ y/3/3tqlo. Thus, we can summarize for the scattering 
from the first fractal iteration: tp < 1 implies entirely coherent 
regime, tp ~ tt and tp ~ 2ir correspond to minimum and max- 
imum, respectively, and finally, the condition tp > 2ir leads to 
entirely incoherent regime. 

The behavior of the structure factor for arbitrary mth iter- 
ation can be described in the same manner. Each of the most 
pronounced minima or maxima corresponds to the interfer- 
ence of cluster amplitudes, where each cluster is a fractal iter- 
ation of the order k = 1, • • • , m. The most common distances 
between the center of masses of the clusters are equal to (see 
Sec.imj 

b k = VzioPt^- 1 . (34) 

Then by analogy with the above considerations, the minima 
and maxima positions can be estimated from the condition that 
the distances (l34l equal n/q and 27r/g, respectively. There- 
fore, we obtain the conditions of minima 

^"Tl^F 1 ' fc = 1 '"-' m ' (35) 

and the conditions of maxima 

^"vW" 1 ' fc = 1 '"-< m - (36) 



These relations are satisfied with a good accuracy, see Fig. [2] 
Thus, the number of minima is equal to the number of fractal 
iteration. The same is true for the number of maxima. 

To explain the generalized power-law behavior, one can 
use the approximate expression d28"T >. which works quite well 
(Fig.Et). Note that for a given range l/(v / 3/3t/3 s fe " 1 ^o) < q < 
l/iVZptPah), only one term SiG8* -1 t?) in the product © 
has a non-trivial behavior. The other terms are nearly equal 
N± = 9 (the coherent regime for the terms on the right) or 
approximately equal one (the incoherent regime for the terms 
on the left). This physically means that the diffraction pattern 
in this range is produced by the interference of only one group 
of subunits, which are the fcth fractal iterations. Hence, by in- 
creasing the argument 1//3 S times, one term on the right in the 
product is forced from the coherent regime into the incoher- 
ent one, thus reducing the total value of the structure factor 
Ni = 9 times. Therefore, we arrive at the equation 

S m (q/l3 s )~p°S m (q), (37) 

valid for 1/(VW ) < q < l/(\/3/3 t /? s m ~%). Here we use 
the equation Ni = l/P^, which follows from the definition 
of fractal dimension (fT9l . Equation d37l l in conjunction with 
the relation d33l l yields 

I m (q/I3 s )(q/P S ) D ~I m (q)q D (38) 

for the same region. This means that within the region, 
the functions S m (q)q D and I m (q)q D are approximately log- 
periodic with the period, equal to the inverse scaling factor 
l//3 s [seeFig.[2>]. 

3. Porod region 

Beyond the fractal region 

2 > lA/WV), (39) 

S m (q) — 1, as discussed above, and we have in accordance 
with Eq. d26]l 

I m (q)/I m (0) - \F (PTqlo/2)\ 2 /N m . (40) 

It follows that beyond the fractal region, the scattering in- 
tensity resembles the intensity of the initiator, i.e., a ball in 
the present case. In p articular, the maxima of the curve obey 
Porod's law q~ 4 I22ll23ll in the Porod region 

q > i/(/3 s m W, (4i) 

whose lower border is of order the upper edge of fractal re- 
gion. Such a behavior can be seen in Fig. [2^. However, if the 
radius of initiator is much smaller than Iq then we have an in- 
termediate region, where I m (q)/Im(0) — 1/N m . So, a well 
or even "shelf appears in the normalized intensity between 
the fractal and Porod regions near the volume 1 /N m j46Tl . 

The generalized power-law behavior of GSSVF is consis- 
tent with the SAS from othe r typ es of deterministic fractals 
reported in the literature J27l H SBH]. In the paper JH], 
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the authors studied SAS from Menger sponges and another 
type of fractals and found some general scattering properties 
of deterministic fractals: the generalized power-law behavior 
of the intensities and log-periodicity of the main maxima and 
minima. However, the method of calculations does not allow 
the authors to distinguish between the fractal structure fac- 
tor and the total scattering intensity; for instance, the correct 
asymptotics of structure factor was not found, and the posi- 
tions of the minima and maxima were not estimated. In the 
paper [43], self-affine characteristics of the scattering curves, 
like Eq. d37l ), were investigated; however, the author studied 
the exact anisotropic formfactors like Eq. d22i >, which were 
not averaged over all directions of momentum. We believe 
that the log-periodicity, even approximate, takes place only 
after averaging the scattering curves over the angles. We dis- 
cuss this important issue in Sec. IVTlbelow. 



V. POLYDISPERSE FORMFACTOR AND STRUCTURE 
FACTOR 

In a real physical system, scatterers almost always have dif- 
ferent sizes. Therefore, a more realistic description should in- 
volve polydispersity . This means that we deal with a set of 
fractals with various sizes and forms, in general. Here we 
can consider a kind of polydispersity, an ensemble of GSSVF 
with different sizes I (that is, I is the length of initial cube, 
see Sec. ITlTb . Note that in the previous sections we denote 
the fractal size as Iq, while here and below Iq is the mean 
value of fractal sizes over the ensemble. The distribution func- 
tion -Dn(0 of the scatterer sizes is defined in such a way that 
Dj$(l)dl gives the probability of finding a fractal whose size 
falls within the interval (1,1 + dl). We consider here one of 
the most common distribution function, the log-normal distri- 
bution, given by 



£>n(0 = 



1 



a/(27r) 1 /2 



exp 



;io g (/Ao) + ( 7 2 /2]- 

2a 2 



, (42) 



where a = [log(l + of )] a . The quantities Iq and cr r are the 
mean length and its relative variance 

,2\l/2 , 



l = (l) D , a T ={{P) D -llY'*/l Q , 



(43) 



where (■■■) D = / °° • • • D N (l)dl. 

As in the previous sections, we assume that spatial positions 
of different fractals are uncorrelated. Hence, the resulted in- 
tensity is the average of the intensity (O over the distribution 
function d42l 

/>oo 

I m (q)=n\Ap\ 2 (\F m (q)\ 2 )V^(l)D N (l)dl, (44) 
Jo 

where the volume and amplitude of monodisperse fractal are 
given by Eqs. ( f2Qb and ( 1221 . respectively. The fractal structure 
factor in the presence of polydispersity can be calculated in 
the same manner but without the term V^(7) 



S m (q) = N r 



r°° i m \ 
. (II \Gi(q)\ 2 )D N mh (45) 

*' i=l 
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FIG. 3. (color online) Influence of polydispersity on the scattering 
curves: (a) intensity d44b and (b) structure factor J45b for various val- 
ues of the relative variance <r r , given by Eq. d43t . Scattering curves 
for cr r = 0.2 and 0.4 are scaled-up for clarity by 10 and 100, respec- 
tively, (c) The quantity I(q)q D clearly shows the log-periodicity in 
the fractal region. The period in the log-scale is equal to log 10 (l//3 s ), 
where /3 S is the scaling factor of fractal. The colors are the same as 
in Fig. [2] 
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where Gi is given by Eq. ( f2TT >. Equations ( f44l > and (145} corre- 
spond to the different statistical ensembles. Indeed, the former 
averaging assumes that the size of balls, constituting the frac- 
tal, equals 1/2 and thus it is also distributed in accordance with 
Eq. d42i >. The latter averaging assumes that all the distances 
between the ball centers inside the fractal are changed when 
the fractal size varies, but the radius of the balls is invariable. 
In this case, the structure factor becomes clearly defined even 
in the presence of polydispersity. 

Figure [3] shows the influence of polydispersity on the ob- 
tained scattering intensity d44t and structure factor d45l ). As 
expected, the small oscillations, present in the monodisperse 
case (see Fig. [2}, are now smeared out, and the scattering 
curves became smoother. One can see that the smoothness 
increases with growing the width of the distribution function, 
which is controlled by the relative variance a T , Nevertheless, 
the exponent of the generalized power-law dependence re- 
mains unchanged and equal D in the fractal region d32l ). In the 
Porod region d4Tb . the value of exponent still equals —4. The 
structure factor is close to one beyond the fractal region d39b . 
The logarithmic periodicity of the quantity I(q)q D becomes 
even more pronounced in comparison with the monodisper- 
sive fractals, discussed in Sec. lIV B~2l 



VI. PAIR DISTRIBUTION FUNCTION AND PAIR 
DISTANCE DISTRIBUTION FUNCTION 

A. General definitions 

The formfactor and related structure factor are intimately 
connected to spatial correlations of subunits, composing the 
investigated system, and the main problem of S AS is to extract 
information about the real-space structure from the scattering 
data. 

Let us consider first a monodisperse fractal, composed of 
balls of the same radius. The general expression ( TT3} can be 
rewritten in the form 



sin gr jk 
qr jk 



(46) 



where rjk = \rj — ru\ are the relative distances between the 
ball centers. Deriving this last relation, we use the formula 
(exp(ia • q)) = sin(aq) / (aq) , which follows from equa- 
tion (O. In Eq. ( l46l . the sum contains many equals terms, 
because the distances separating different points can coincide. 
One can introduce the probability density of finding the dis- 
tance r between the centers of two arbitrarily taken balls in- 
side the fractal: 



p(r) 



r - r jk j 



N(N - 



— ^C p 5{r -r p ), 



where C p are the numbers of distances separated by r p . 
dependence on the fractal length I appears through r p = 



(47) 

The 

I z ri . 




FIG. 4. (Color online) (a) The coefficients C v in the expression 
for the pair distribution function of GSSVF. The centers of groups 
are the distances between the clusters d34t . and the period on the log- 
scale is equal to log 10 (l//3 s ). (b) The influence of polydispersity on 
the pair distance distribution function. Scattering curves for the dis- 
tribution width <7r = 0.2, 0.4 are scaled-up for clarity by 10 and 10 2 
respectively. Solid (black) lines show the solutions obtained with 
Eq. {49} from the polydisperse structure factor of Fig. [3p, dashed 
(blue) lines represent its direct calculation {52}. The shift of max- 
ima in comparison to Eq. {34} is due to a polydispersity effect. The 
exponent of the generalized power-law behavior for p(r) is equal to 
D — 1, see Eq. {54}- ( c ) The log-periodic function N(r)/r D with 
the same log-period. The quantity N(r) {50} is proportional to the 
"mass", or measure, enclosed in the imaginary sphere of radius r, 
which is centered on fractal. 
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where the distances z p correspond to the fractal of unit length. 
The quantity ( |47t is called pair distance distribution function. 
In terms of this function, Eq. d4oT ) reads 



S(q) = 1 + (N-1) 



+00 



drp(r) 



sin qr 
qr 



(48) 



By performing the inverse Fourier sine transform, we obtain 



p{r) 



2 

7T ./II 



S(q) ~ 1 

N — 1 



qr sin qr dq. 



(49) 



Let us choose an arbitrary ball in the fractal and consider 
a spherical layer of radius r and width dr, whose center co- 
incides with the center of the chosen ball. It follows from 
the definition (|47| | that dN = (N — l)p(r)dr gives the aver- 
age number of other ball centers, which lay within the spheri- 
cal layer. Then the total average number of other balls in the 
sphere of radius r is given by 



N(r) = (N - 1) / p(r')dr'. 
Jo 



(50) 



Conversely, (N — l)p(r) = dN(r)/dr. The quantity N(r) 
( TSUI is proportional to the "mass", or measure, enclosed in the 
imaginary sphere of radius r, which is centered on fractal, see 
the discussion in Appendix[A] The relation ( T50b gives a more 
precise definition of N(r), discussed in SecU By substituting 
Eq. d49l into Eq. d50l ), we arrive at the relation between N(r) 
and the fractal structure factor 

2 f +oc S(q) — 1 
N(r) = — / [sin(qr) — qr cos(qr)] dq. (51) 



The polydispersity can be taken into account by averaging 
the above equations d47li-(l5l1l over the distribution (l42l . All 
the relations take the same form but with the replacements 

S(q) -► (S(q)) D , N(r) -> (N(r)) D ,p(r) -> (p(r)) D . The 
last quantity can be calculated explicitly 



(p(r)) D = 



1 ) Zr> \ Z n 



N(N - 1) 



(52) 



For a finite iteration, the total number of balls N is given by 
Eq. (H7j. 

For describing spatial particle correlations, one can also use 
pair distribution function pill g(r), which is directly related 
to the pair distance distribution function d47l ) 



_ p(r)l 3 



4-71-7.2 



(53) 



where Z 3 is the total volume of fractal. As discussed in Sec. [H 
the pair distribution function has a transparent physical in- 
terpretation: it is nothing else but the conditional probabil- 
ity density, because it gives the probability density to find a 
particle at the distance r from another particle, provided a po- 
sition of the latter particle is given. As follows from the dis- 
cussion below Eq. (1491 , the number of particles within the 
spherical layer is indeed given by dN = ng(r)dV. Here 



n = (N — I) /I 3 is the average particle density in the frac- 
tal (if we neglect the difference between N and N — 1 for 
a large number of particles), and dV^ = 47rr 2 dr is the layer 
volume. 

Note that Eq. (1531 is consistent with the other definition 
g(r) = l 3 /[N(N — 1)1 (y^ y S( r — r^)), more common in 
the literature (see, e.g., Ref. 14711 ). Here the brackets denote 
the average over all the directions of r, and ry are the rel- 
ative positions of two particles. By using the average of the 
(5-function (6(r — a)) = 5(r — a)/(4irr 2 ), we derive Eq. d53l>. 



B. Analysis of results 

The real-space characteristics are shown in Fig. |4] Figure 

represents the coefficients in the expression ( f4Tb for the pair 
distance distribution function. They are found numerically for 
the fourth iteration by a simple combinatoric analysis. The 
self-similarity of GSSVF manifests itself in the periodicity of 
the groups of distances on the logarithmic scale. As expected, 
the centers of groups are the distances between the clusters 
(l34i >. discussed in Sec. lIVB~2l and the period in the log-scale 
is related to the scaling parameter through log 10 (l/ /3 S ). 

When the polydispersity is present, the pair distance distri- 
bution function can be calculated from Eq. d49l with the poly- 
disperse structure factor (l45l l or directly from Eq. d52l with 
known quantities C p and z p . We obtain a very good agreement 
between the two formulas, see Fig.|4j). Small differences ap- 
pear at the low distances because of the finite upper integration 
limit in Eq. ( 149) , which we have to use in numerical calcula- 
tions. Equation ( |49l is nevertheless preferable for large fractal 
iterations, since it surpasses the inconvenience of calculating 
the exponentially increasing number of distances. This gives 
a nearly independent computation time with increasing the it- 
eration number. 

The simple power law (0]i for g(r) implies that p(r) oc 
r ^ 1 due to the relation ( [S3) . Instead, one can see the gen- 
eralized power-law behavior in Fig. [4}}. Moreover, the quan- 
tity p(r) I r D ~ x is log-periodic in the fractal region in the real 
space, which results from the self-affinity relation ( |37| ) for the 
structure factor. Indeed, the main contribution in the integral 
of Eq. d49l comes from the fractal region in momentum space, 
where S(q) 1 and Eq. ( f37b is satisfied. Then, equating 
S(q) — 1 to S(q) in the integral of Eq. ( |49l and substituting 
q^ q/(3 s yield 



p(/3 s r) p(r) 



(M 1 



(54) 



which is valid for /3 t ^- 2 l < r < l . By Eq. (|53}, the anal- 
ogous relation takes place for the pair distribution function 



g(j3 s r) g(r) 



D-3 r D-3 ' 



(55) 



By using Eq. ( BTT l in the same manner as Eq. d49l ), we arrive 
at the relation 



N((3 s r) = N(r) 
(f3 s r) D r D ' 



(56) 
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The last three equations are equivalent. 

The real-space correlation functions of fractal measure are 
well-studied in the literature, see the textbooks lfl2l l33tl and 
references therein. In particular, the log-periodicity of the re- 
duced mass-radius relation d56l ) is known for ideal determin- 
istic fractals (see, e.g., Refs. II 52145411 . In this paper, the mass- 
radius relation is connected to the log-periodicity ( f38T > in the 
momentum space for finite fractal iterations. Conversely, once 
Eq. ( 154-b or 05] ) or (T5~6b is fulfilled, then the log-periodicity 
S(q)q D , given by Eq. d37l ). is also satisfied. For a fractal em- 
bedded into the 3D space, the quantity N(r) is radially sym- 
metric by definition. This means that the log-periodicity of 
S(q)q D , which follows from Eqs. (f5TT > and d56l ), is met only 
for the structure factor (|25l l, which includes the averaging over 
all directions of momentum. The squared anisotropic scat- 
tering amplitude YiiLi \Gi(q)\ 2 does not obey in general the 
self-affinity relation (|571 l, as discussed in Sec. II V B| 

The log-periodicity of the mass-radius relation d56l > is 
shown in Fig. 2J:. The number of periods on the logarithmic 
scale coincides with the fractal iteration number. 

Polydispersity makes the above log-periodicity even more 
apparent. Increasing the value of the distribution width a r 
leads to smoothing of p(r), g(r) and N(r). One can see a 
little shift of the minima to the right. From the general con- 
siderations, this is allowed as long as the distribution width is 
much greater than the relative shift Ar/r. The simple power 
law [p(r) oc r ^ 1 , g(r) oc r D ~ 3 , and N(r) oc r D ], typical 
for random fractals, is restored in the limit of strongly devel- 
oped polydispersity. 



VII. CONCLUSION 



general features, common for deterministic mass fractals with 
a single scale, and can be used for interpreting the experimen- 
tal data. For this kind of fractal, one can extract a number of 
parameters from the scattering intensity, see Fig. [3] 

( i) The fractal dimension from the generalized power law. 

( ii) The fractal scaling parameter from the period on the loga- 
rithmic scale. 

( Hi) The number of fractal iteration, which is equal to the num- 
ber of periods of function I(q)q D . 

(iv) The lower and upper fractal edges from the diagram 
I(q)q D as the beginning and end of "periodicity region". 
They allow us to estimate the fractal size and the smallest dis- 
tance between fractal units in accordance with Eq. (01. 

(v) The total number of structural units, from which the fractal 
is composed, by the relation N m — (1/ (3 s ) mD . 

Note that the number of structural units can also be estimated 
from the ratio /(0)//(g max ), where q max is the upper edge. 

Once the internal fractal structure is known in more detail, 
one can derive and use an analytical expression for the scat- 
tering amplitude such as Eq. (l22l and use it directly to fit the 
scattering data. This scheme could be applied even more ef- 
ficiently to anisotropic scattering with a position sensitive de- 
tector. 

The results obtained can be applied for various structures, 
whose geometries are based on iterations of fractal systems. 
This includes magnetic cluster structures, artificially created 
chemical compounds and so on. 

We consider mass fractals composed of the same units, but 
the developed scheme allows a generalization to mass fractals 
containing units of various shapes and sizes and to surface 
fractals as well. Generalizations of the developed scheme are 
under way. 



We develop a model of deterministic fractal that enables 
explicit analytic solutions for the scattering amplitude. The 
model generalizes the regular 3D Vicsek fractals. The proper- 
ties of GSSVF are studied and analyzed in the both momen- 
tum and real spaces. The fractal dimension is controlled by the 
scaling parameter and can vary from to 2.862 . . .. Based on 
the developed model, we derive analytical expressions for the 
main properties in the monodisperse and polydisperse cases: 
the formfactor and fractal structure factor, the pair distribution 
function, the pair distance distribution function, radius of gy- 
ration, intensity in zero angle and the edges of fractal region 
in the momentum and real spaces. We found the logarithmic 
periodicity of the intensity scaled q D times within the fractal 
region [see Eq. d38l >l and relate it to the log-periodicity of the 
functions (l54l i, (1551 . and (l56l l within the fractal region of real 
space. The period is governed by the scaling parameter (3 S 
and equal to log 10 (l//3 s ) on the logarithmic scale. This be- 
havior of the scattering curves is explained and interpreted by 
analogy with optics (see Sec. lIVB~2l . 

An advantage of the exact analytic solutions for determin- 
istic fractals is the absence of phenomenological parameters, 
often used in the random fractal models. The behavior of ran- 
dom fractals can be simulated by introducing strong polydis- 
persity in the deterministic fractal models. 

The results, obtained for the GSSVF, illustrate a number of 
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Appendix A: Hausdorff dimension 

The Hausdorff dimension lf55]l can be rigorously defined as 
follows (see, e.g., Refs. Ill2il33l 15611 ). 

Let A be a subset of the n-dimensional Euclidean space and 
{Vi} be a covering of A with = diam(T^) ^ a. Then the 
a-dimensional Hausdorff measure m a (A) of the set A is 

m a (A) ee lim inf V a? , a> 0. (Al) 

a^0 {V z } ^ 
i 

Here the infimum is on all possible covering. Note that 
m a (A) can be infinite and a is not integer in general. 
We define the Hausdorff dimension D of the set A by 

D ee inf{a : m a {A) = 0} = sup{a : m a (A) = +oo}. 

(A2) 
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In other words, the Hausdorff dimension is the value of a for 
which the Hausdorff measure jumps from zero to infinity. For 
the value a = D, this measure may be anywhere between 
zero and infinity. 

In practice, the rigorous definition iA2\ is rather difficult to 
apply, and one can use other methods for calculating the Haus- 
dorff (fractal) dimension of a fractal 11121 l33fl . For instance, 
one can use the mass-radius relation, that is, the "mass" of the 
structure within a ball of dimension n and radius r centered 
on the fractal 

M(r) = A(r)r D , (A3) 

where logA(r)/logr — > as r —> oo. By mass we mean 
the total fractal measure, which could be a mass, volume, sur- 
face area or any other scalar quantity attached to the fractal 
support. 

Let us consider a few examples. For a self-similar deter- 
ministic fractal of total length L, whose first iteration consists 
of k elements of size f3 s L, one can write M(L) = kM((3 s L). 
Using Eq. ( lA31 l. we obtain 

fc/3f - 1. (A4) 

The Vicsek fractal, considered in Sec. [Ell] is a particular case 
of the above fractal, and the formula (fT9l i for its fractal dimen- 
sion results from Eq. ( 1A41 > at k = 9. For a multiscale fractal, 
giving at each iteration fcj elements of size (3 S iL, we obtain 

X>/3£ = L. (A5) 



A more complicated example is a dense circle packing. It 
is an infinite set of non-overlapping circles of smaller and 
smaller radii inscribed into a larger circle in order to fill it 
completely. The set of circles is obviously a fractal. The dis- 
tribution n(r) of the circle radii is given by a simple power 
law n(r) oc r~ T with 2 < r < 3 (see, e.g., Ref. [15711 ). The 
exponent r can easily be related to the fractal dimension. An 
analog of finite iteration is the cutoff length a, for which only 
the circles of radii larger than a are considered. The minimal 
number of disks of radius a needed to cover a circle of radius 
r is proportional to r/a. Then the minimal number of disks 
for covering the fractal with a finite cutoff length a is 

p oo 

N(a) oc / drn(r)r/a oc a 1_r . (A6) 

J a 

Comparing this equation with the definition of Hausdorff di- 
mension N(a) oc a~ D yields D = r — 1. The value of dimen- 
sion depends on a specific type of packing and can be found 
numerically. For instance, the value D — 1.307. . . is ob- 
tained for the classical Apollonian packing 15711 . Note that if 
we consider the set of the filled circles (disks), then its Haus- 
dorff dimension is obviously equal to D = 2. Indeed, the 
minimal number of disks of radius a needed to cover a disk of 
radius r is proportional to r 2 /a 2 , and we obtain in the same 
manner 

/>oo 

N(a) oc / drn(r)r 2 /a 2 oc a~ 2 , (A7) 

J a 

because the integral of n(r)r 2 converges when a —> for 
r < 3. 
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